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There is very strong circumstantial evidence that there was an inflationary 
epoch very early in the history of the universe. In this lecture I will describe 
how we might be able to piece together some understanding of the dynamics 
during and immediately after the inflationary epoch. 



(N, 
> ■ 1 Introduction 

cn ' We live in a very large, very old, and nearly (perhaps exactly) spatially fiat universe. The 

^ ■ universe we observe is, at least on large scales, remarkably homogeneous and isotropic. 

Q^ . These attributes of our universe can be ascribed to a period of very rapid expansion, or 

^ ! inflation, at some time during the very early history of the universe |I| . 

1-^ I A sufficiently long epoch of primordial inflation leads to a homogeneous/isotropic 

I I universe that is old and flat. The really good news is that very many reasonable models 

tj'' have been proposed for inflation.^ In some ways, inflation is generic. That is also the 

^ . really bad news, since we would like to use the early universe to learn something about 

physics at energy scales we can't produce in terrestrial laboratories. We want to know 
'k/< I more than just that inflation occurred, we want to learn something of the dynamics of 

^ I the expansion of the universe during inflation. That may tell us something about physics 

at very high energies. It would also allow us to restrict the number of inflation models. 
If we can differentiate between various inflation models, then inflation can be used as a 
phenomenological guide for understanding physics at very high energies. 

Oldness, flatness, and homogeneity/isotropy only tell us the minimum length of the 
inflationary era. They are not very useful instruments to probe the dynamics of inflation. 
If that is our goal, we must flnd another tool. In this lecture I will discuss two things 
associated with inflation that allow us to probe the dynamics of inflation: perturbations 
and preheating/reheating. 

While the universe is homogeneous and isotropic on large scales, it is inhomogeneous 
on small scales. The inhomogeneity in the distribution of galaxies, clusters, and other 

^Perhaps a more accurate statement is that there are many models that seemed reasonable to the 
people who proposed them at the time they were proposed. 



luminous objects is believed to result from small seed primordial perturbations in the 
density field produced during infiation. Density perturbations produced in infiation also 
lead to the observed anisotropy in the temperature of the cosmic microwave background 
radiation. A background of primordial gravitational waves is also produced during infia- 
tion. While the background gravitational waves do not provide the seeds or infiuence the 
development of structure, gravitational waves do lead to temperature variations in the 
cosmic microwave background radiation. If we can extract the primordial density per- 
turbations from observations of large-scale structure and cosmic microwave background 
radiation temperature fiuctuations, then we can learn something about the dynamics of 
infiation. If we can discover evidence of a gravitational wave background, then we will 
know even more about the dynamics of the expansion rate during infiation. 

Infiation was wonderful, but all good things must end. The early universe somehow 
made the transition from an infiationary phase to a radiation-dominated phase. Perhaps 
there was a brief matter-dominated phase between the infiation and radiation eras. In the 
past few years we have come to appreciate that some interesting phenomena like phase 
transitions, baryogenesis, and dark matter production can occur at the end of infiation. 
Perhaps by studying how infiation ended, we can learn something of the dynamics of the 
universe during infiation. 

2 Perturbations Produced During Inflation 

One of the striking features of the cosmic background radiation (GBR) temperature 
fiuctuations is the growing evidence that the fiuctuations are acausal.0 The GBR fiuc- 
tuations were largely imprinted at the time of last-scattering, about 300,000 years after 
the bang. However, there seems to be fiuctuations on length scales much larger than 
300,000 light years!^ How could a causal process imprint correlations on scales larger 
than the light-travel distance since the time of the bang? The answer is infiation. 

In order to see how infiation solves this problem, first consider the evolution of the 
Hubble radius with the scale factor a(t):0 
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"' < a" (RD) 

x3/2 (MD). 

In a fc = matter-dominated universe the age is related to if by t = {2/3)H~^, so 
Rh = (3/2)t. In the early radiation-dominated universe t = {1/2)H~^, so Rh = 2t. 

On length scales smaller than Rh it is possible to move material around and make an 
imprint upon the universe. Scales larger than Rh are "beyond the Hubble radius," and 

^Exactly what is meant by acausal will be explained shortly. Acausal is in fact somewhat of a 
misnomer since, as we shall see, inflation produces "acausal" perturbations by completely causal physics. 

^Although definitive data is not yet in hand, the issue of the existence of acausal perturbations will 
be settled very soon. 

''Here and throughout the paper "RD" is short for radiation dominated, and "MD" implies matter 
dominated. 
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Figure 1: Physical sizes increase as a{t) in the expanding universe. The Hubble radius 
evolves as Rh = H~^ = (87rGp(a)/3)^/^. In a radiation-dominated or matter-dominated 
universe (illustrated by the left panel) any physical length scale A starts larger than R^, 
then crosses the Hubble radius (A = H~^) only once. However, if there was a period of 
early inflation when Rh increased more slowly than a (as illustrated in the right panel), 
it is possible for a physical length scale to start smaller than Rh, become larger than 
Rh, and after inflation ends become once again smaller than Rh- Periods during which 
the scale is larger than the Hubble radius are indicated by the dotted line. 



the expansion of the universe prevents the establishment of any perturbation on scales 
larger than Rh- 

Next consider the evolution of some physical length scale A. Clearly, any physical 
length scale changes in expansion in proportion to a{t). 

Correlations on physical length scales A larger than Rnit) are often called acausal. 

Now let us form the dimensionless ratio L = X/Rh- If L is smaller than unity, 
the length scale is smaller than the Hubble radius and it is possible to imagine some 
microphysical process establishing perturbations on that scale, while if L is larger than 
unity, no microphysical process can account for perturbations on that scale. 

Since Rh = a/ a, and A oc a, the ratio L is proportional to a, and L scales as a, which 
in turn is proportional to — (p + 3p). There are two possible scenarios for L depending 



upon the sign of p + 3p: 

• J < — > Rh grows faster than A, happens for p + 3p > , . 

I > ^ Rh grows more slowly than A, happens for p + 3p < 0. 

In the standard MD or RD universe, p + 3p > 0, and Rh grows faster than A. 

Perturbations that appear to be noncausal at last scattering can be produced if 
sometime during the early evolution of the universe the expansion was such that a > 0. If 
a > 0, then the Hubble radius will increase more slowly than any physical scale increases 
during inflation. Physical length scales larger than the Hubble radius at the time of last 
scattering would have been smaller than the Hubble radius during the accelerated era. 
Therefore it is possible to imprint correlations during inflation as a scale passes out of 
the horizon and have it appear as an acausal perturbation at last scattering. 

From Einstein's equations, the "acceleration" is related to the energy density and the 
pressure as a oc — (p + 3p). Therefore, an acceleration (positive a) requires an unusual 
equation of state with p + 3p < 0. This is the condition for "accelerated expansion" or 
"inflation." 

We do not yet know when inflation occurred, but the best guess for the different 
epochs in the history of the universe is given in Table |l|. It is useful to spend a few 
minutes discussing the movements of Table |I|. 

The first movement of the Cosmic Symphony may be dominated by the string section 
if on the smallest scales there is a fundamental stringiness to elementary particles. If this 
is true, then the first movement in the cosmic symphony would have been a pizzicato 
movement of vibrating strings about 10~'^^s after the bang. There is basically nothing 
known about the stringy phase, if indeed there was one. We do not yet know enough 
about this era to predict relics, or even the equation of state. 

The earliest phase we have information about is the inflationary phase. The infla- 
tionary movement probably followed the string movement, lasting approximately 10~^^ 
seconds. During inflation the energy density of the universe was dominated by vacuum 
energy, with equation of state pv — —pv- As we shall see, the best information we have 
of the inflationary phase is from the quantum fluctuations during inflation, which were 
imprinted upon the metric, and can be observed as CBR fluctuations and the departures 
from homogeneity and isotropy in the matter distribution, e.g., the power spectrum. 
Inflation also produces a background of gravitational radiation, which can be detected 
by its effect on the CBR, or if inflation was sufficiently exotic, by direct detection of the 
relic background by experiments such as LIGO or LISA. 

Inflation was wonderful, but all good things must end. A lot of effort has gone into 
studying the end of inflation (For a review, see Kofman et al., 0.) It was likely that 
there was a brief period during which the energy density of the universe was dominated 
by coherent oscillations of the inflaton field. During coherent oscillations the infiaton 
energy density scales as a~^ where a is the scale factor, so the expansion rate of the 
universe decreased as in a matter-dominated universe with p ~ 0. Very little is known 
about this period immediately after infiation, but there is hope that one day we will 
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Table 1: Different epochs in the history of the universe and the associated tempos of 
the ever decreasing expansion rate H, along with the equation of state and some of the 
relics produced during the various eras. 



discover a relic. Noteworthy events that might have occurred during this phase include 
baryogenesis, phase transitions, and generation of dark matter. 

We do know that the universe was radiation dominated for almost all of the first 
10,000 years. The best preserved relics of the radiation-dominated era are the light 
elements. The light elements were produced in the radiation-dominated universe one 
second to three minutes after the bang.Q If the baryon asymmetry is associated with the 
electroweak transition, then the asymmetry was generated in the radiation-dominated 
era. The radiation era is also a likely source of dark matter such as WIMPS or axions. 
If one day we can detect the 1.9 K neutrino background, it would be a direct relic of the 
radiation era. The equation of state during the radiation era is pr = Pr/3. 

The earliest picture of the matter-dominated era is the CBR. Recombination and 
matter-radiation decoupling occurred while the universe was matter dominated. Struc- 
ture developed from small primordial seeds during the matter-dominated era. The pres- 
sure is negligible during the matter-dominated era. 

Finally, if recent determinations of the Hubble diagram from observations of distant 
high-redshift Type-I supernovae are correctly interpreted, the expansion of the universe 
is increasing today (a > 0). This would mean that the universe has recently embarked 
on another inflationary era, but with the Hubble expansion rate much less than the rate 
during the first inflationary era. 

2.1 Simple Dynamics of Inflation: The Inflaton 

In building inflation models it is necessary to find a mechanism by which a universe 
dominated by vacuum energy can make a transition from the inflationary universe to a 
matter-dominated or radiation-dominated universe. There is some unknown dynamics 
causing the expansion rate to change with time. There may be several degrees of freedom 
involved in determining the expansion rate during inflation, but the simplest assumption 
is that there is only one dynamical degree of freedom responsible for the evolution of the 
expansion rate. 

If there is a single degree of freedom at work during inflation, then the evolution from 
the inflationary phase may be modeled by the action of a scalar fleld (p evolving under 
the influence of a potential V{(j)). Let's imagine the scalar fleld is displaced from the 
minimum of its potential as illustrated in Fig. ^. If the energy density of the universe 
is dominated by the potential energy of the scalar fleld 0, known as the inflaton, then 
p + 3p will be negative. The vacuum energy disappears when the scalar fleld evolves to its 
minimum. The amount of time required for the scalar fleld to evolve to its minimum and 
inflation to end (or even more useful, the number of e-folds of growth of the scale factor) 
can be found by solving the classical field equation for the evolution of the infiaton field. 



^Although I may speak of time after the bang, I will not address the issue of whether the universe 
had a beginning or not, which in the modern context is the question of whether inflation is eternal. 
For the purpose of this discussion, time zero of the bang can be taken as some time before the end of 
inflation in the region of the universe we observe. 



which is simply given by + 3iJ0 + dV/dip = 0. 

2.2 Quantum fluctuations 

In addition to the classical motion of the inflaton field, during infiation there are quantum 
fiuctuations.0 Since the total energy density of the universe is dominated by the infiaton 
potential energy density, fiuctuations in the inflaton field lead to fiuctuations in the 
energy density. Because of the rapid expansion of the universe during infiation, these 
fiuctuations in the energy density are frozen into super-Hubble-radius-size perturbations. 
Later, in the radiation or matter-dominated era they will come within the Hubble radius 
as if they were noncausal perturbations. 

The spectrum and amplitude of perturbations depend upon the nature of the infiaton 
potential. Mukhanov ||^ has developed a very nice formalism for the calculation of density 
perturbations. One starts with the action for gravity (the Einstein-Hilbert action) plus 
a minimally-coupled scalar infiaton field 0: 



S = — d'^x y/^ 



^i? - lg^''d,<pd,<P + V{<p) 



(3) 



Here R is the Ricci curvature scalar. Quantum fiuctuations result in perturbations in 
the metric tensor and the infiaton field 



9,,^ ^ gjlu + S9,u ; 0^00 + 50, (4) 

where gj^^^ is the Friedmann-Robertson- Walker metric, and <f)o{t) is the classical solu- 
tion for the homogeneous, isotropic evolution of the infiaton. The action describing the 
dynamics of the small perturbations can be written as 



'^25' = - I d'^x 



dr'^ 



aHH , (5) 



i.e., the action in conformal time r [dr"^ = a^{t)dt^) for a scalar field in Minkowski 
space, with mass-squared m^ = —z'^d'^z/dr'^. Here, the scalar field m is a combination 
of metric fiuctuations dg^y and scalar field fiuctuations 54>. This scalar field is related to 
the amplitude of the density perturbation. 

The simple matter of calculating the perturbation spectrum for a noninteracting 
scalar field in Minkowski space will give the amplitude and spectrum of the density 
perturbations. The problem is that the solution to the field equations depends upon 
the background field evolution through the dependence of the mass of the field upon 
z. Different choices for the inflaton potential V^(0) results in different background fleld 
evolutions, and hence, different spectra and amplitudes for the density perturbations. 

Before proceeding, now is a useful time to remark that in addition to scalar density 
perturbations, there are also fluctuations in the transverse, traceless component of the 

^Here I will continue to assume there is only one dynamical degree of freedom. 



spatial part of the metric. These fluctuations (known as tensor fluctuations) can be 
thought of as a background of gravitons. 

Although the scalar and tensor spectra depend upon V"(0), for most potentials they 
can be characterized by Q^ms (^^^ amplitude of the scalar and tensor spectra on large 
length scales added in quadrature), n (the scalar spectral index describing the best 
power-law fit of the primordial scalar spectrum), r (the ratio of the tensor-to-scalar 
contribution to C2 in the angular power spectrum), and tit ( the tensor spectral index 
describing the best power-law fit of the primordial tensor spectrum). For single- field, 
slow-roll inflation models, there is a relationship between tit and r, so in fact there are 
only three independent variables. Furthermore, the amplitude of the fluctuations often 
depends upon a free parameter in the potential, and the spectra are normalized by Qrms- 
This leads to a characterization of a wide-range of inflaton potentials in terms of two 
numbers, n and r. 

2.3 Models of inflation 

A quick perusal of the literature will reveal many models of inflation. Some of the 
familiar names to be found include: old, new, pre-owned, chaotic, quixotic, ergodic, 
exotic, heterotic, autoerotic, natural, supernatural, au natural, power-law, powerless, 
power-mad, one-field, two-field, home-field, modulus, modulo, moduli, self-reproducing, 
self-promoting, hybrid, low-bred, white-bread, first-order, second-order, new-world order, 
pre-big-bang, no-big-bang, post-big-bang, D-term, F-term, winter-term, super symmetric, 
superstring, superstitious, extended, hyperextended, overextended, D-brane, p-brane. 
No-brain, dilaton, dilettante, . . . 

Probably the first step in sorting out different models is a classification scheme. One 
proposed classification scheme has two main types. Type-I inflation models are models 
based on a single inflaton field, slowly rolling under the influence of an inflaton potential 
V^(0). This may seem like a restrictive class, but in fact many more complicated models 
can be expressed in terms of an equivalent Type-I model. For instance "extended" 
inflation, which is a Jordan-Brans-Dicke model with an inflaton field and a JBD scalar 
field, can be recast as an effective Type-I model. Anything that is not Type I is denoted 
as a Type-II model. 

There are subclasses within Type I. Type-la models are "large-field" models, where 
the inflaton field starts large and evolves toward its minimum. Examples of large-field 
models are chaotic inflation and power-law inflation. Type-lb models are "small-field" 
models, where the inflaton field starts small and evolves to its minimum at larger values of 
the inflaton field. Examples of small-field models are new inflation and natural inflation. 
Finally, hybrid-inflation models are classified as Type-Ic models. In Type-la and Type-lb 
models the vacuum energy is approximately zero at the end of inflation. Hybrid models 
have a significant vacuum energy at the end of inflation (see Fig ^. Hybrid inflation is 
usually terminated by a first-order phase transition or by the action of a second scalar 
field. A more accurate description of large-field and small- field potential is the sign of 
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Figure 2: Schematic illustrations of the inflaton potential energy. The first potential 
is a "large- field" model (Type-la). The second figure illustrates a "small-field" model 
(Type- lb). The final figure illustrates an example of hybrid inflation (Type-Ic). Notice 
that the minimum of the potential is nonzero in the hybrid model. 

the second derivative of the potential: large-field models have V" > while small-field 
models have V" < 0. 

Of course a classification scheme is only reasonable if there are some observable 
quantities that can differentiate between different schemes. It turns out that the different 
Type-I models fill in different regions of the n-r plane, as shown in Fig. |^ (from |^). 
For a given spectral index n, small-field models have a smaller value of r. Shown as an 
ellipse is a very conservative estimate of the uncertainties in n and r that are expected 
after the next round of satellite observations. Although we don't know where the error 
ellipse will fall on the graph, an error ellipse at the indicated size will restrict models. 
So we expect a true inflation phenomenology, where models of inflation are confronted 
by precision observations. 



2.4 Inflation Models in the Era of Precision Cosmology 

It was once said that the words "precision" and "cosmology" could not both appear in a 
sentence containing an even number of negatives. However that statement is now out of 
date, or at the very least, very soon will be out of date. A number of new instruments 
will come on line in the next few years and revolutionize cosmology. 

There is now a world-wide campaign to pin down the microwave anisotropics. In 
the near future, long-duration balloon flights, as well as observations from the cold, 
dry observatory in Antarctica will completely change the situation. Finally, in the next 
decade two satellites, a NASA mission — the Microwave Anisotropy Probe (MAP) — 




Figure 3: Large-field (Type-la), small- field (Type-lb), and hybrid (Type-Ic) models 
seem to populate different regions of the n-r plane. Observational determination of n 
and r could rule out entire classes of models. The dark ellipse indicated the anticipated 
size of the errors in the post-Planck era. Of course the location is yet to be determined. 
This figure is from Dodelson, Kinney and Kolb Q. 

and an ESA mission — PLANCK, will culminate in a determination of the spectrum with 
errors much smaller than the present errors. Of course we don't know what the shape of 
the spectrum will turn out to be, but we can anticipate errors as small as shown in the 
figure. 

With errors of this magnitude, fitting the spectrum will allow determination of n and 
r, useful for inflation, as well as determination of the cosmological parameters {Hq, VLq, 
Qb, a, etc.) to a few percent. 

2.5 Reconstruction 

In addition to restricting the class of inflation models (Type la. Type lb, etc.), it may be 
possible to use the data from precision microwave experiments to reconstruct a fragment 
of the inflationary potential. 

Reconstruction of the inflaton potential (see Lidsey, et al., |^ for a review) refers to 
the process of using observational data, especially microwave background anisotropies, to 
determine the inflaton potential capable of generating the perturbation spectra inferred 
from observations . Of course there is no way to prove that the reconstructed inflaton 
potential was the agent responsible for generating the perturbations. What can be hoped 
for is that one can determine a unique (within observational errors) inflaton potential 
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Table 2: A summary of the slow-roll parameters, H, e, ry, and ^, needed to reconstruct 
a given derivative of the potential to a certain order. Note that the next-order result is 
exact. 

capable of producing the observed perturbation spectra. The reconstructed inflaton 
potential may well be the first concrete piece of information to be obtained about physics 
at scales close to the Planck scale. 

As is well known, inflation produces both scalar and tensor perturbations, and each 
generate microwave anisotropics (see Liddle and Lyth |^ for a review). If V{(f)) is known, 
the perturbation spectra can be computed exactly in linear perturbation theory through 
integration of the relevant mode equations [|]. If the scalar field is rolling sufficiently 
slowly, the solutions to the mode equations may be approximated using something known 
as the slow-roll expansion P, |10|, |ll| . The standard reconstruction program makes use of 
the slow-roll expansion, taking advantage of a calculation of the perturbation spectra by 

m 



which gives the next-order correction to the usual lowest-order 



Stewart and Lyth 
slow-roll results. 

Two crucial ingredients for reconstruction are the primordial scalar and tensor per- 
turbation spectra As{k) and AT{k) which are defined as in Lidsey et ai, p. The scalar 
and tensor perturbations depend upon the behavior of the expansion rate during infia- 
tion, which in turn depends on the value of the scaler infiaton field 0. In order to track 
the change in the expansion rate, we define slow-roll parameters e, r], and ^ as 
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where here the prime superscript implies d/dcf). If the slow-roll parameters are small, e, 
7], and ^ can be expressed in terms of derivatives of the infiaton potential: 
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Table 3: The inflation parameters may be expressed in terms of the primordial scalar 
and tensor perturbation spectra y4|, A^, and the scalar and tensor spectral indices n, 
and dn/d Ink. 

As long as the slow-roll parameters are small compared to unity, the scalar and tensor 
perturbation amplitudes As{k) and Axik) are given by (see Stewart and Lyth [|l2l for 
the normalization) 
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where H, e, and i] are to be determined at the value of when k = aH during inflation, 
and where C = — 2-|-ln2-|-7 ~ —0.73 is a numerical constant, 7 being the Euler constant. 
These equations are the basis of the reconstruction process. 

From the expressions for A'g and A"^ one can express the scalar and tensor spectral 
indices, defined as 
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in terms of slow-roll parameters. 

Perturbative reconstruction requires that one fits an expansion, usually a Taylor series 
of the form 



k 
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where n* = n{k^,), to the observed spectrum in order to extract the coefficients, where 
stars indicate the value at k^,. The scale k^, is most wisely chosen to be at the (logarithmic) 
center of the data, about A;* = 0.01 h Mpc~^. 

For reconstruction, one takes a Hamilton-Jacobi approach where the expansion rate 
is considered fundamental, and the expansion rate is parameterized by a value of the 
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Table 4: A summary of the observables needed to reconstruct a given derivative of the 
potential to a certain order. R^ is the ratio of the tensor to scalar contributions to the 
CBR anisotropics at any conveniently chosen scale. Qrms i^ ^^^ normalization of the 
total (scalar plus tensor) contributions at small /, n* and n'^ refer to the scalar spectral 
index and its derivative at some scale. Knowledge of R^ and Q^ms ^^ equivalent to 
determination of A^lk^) and Aglk^). 



inflaton field. The Friedmann equation may be expressed in terms of the potential V{(f)), 
the expansion rate H{(f)), and the derivative of H as 
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or equivalently in terms of the slow-roll parameter e, as 
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Subsequent derivatives of V may be expressed in terms of additional slow-roll parameters: 
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and so on. So if one can determine the slow-roll parameters, one has information about 
the potential. The slow-roll parameters needed to reconstruct a given derivative of the 
potential is given in Table ^ 

Of course, the problem is that the slow-roll parameters are not directly observable! 
But one can construct an iterative scheme to express the slow-roll parameters in terms 
of observables. The result is shown in Table |^. 

Using the notation ra'^ to indicate d7i/dlnk\^ and i?* = A'^{k^)/A'g{k^:), the recon- 
struction equations are 
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The observables needed to reconstruct a given derivation of the potential are listed in 
Table |. 

The biggest hurdle for successful reconstruction is that many inflation models pre- 
dict a tensor perturbation amphtude (hence, R^) well below the expected threshold for 
detection, and even above detection threshold the errors can be considerable. If the 
tensor modes cannot be identifled a unique reconstruction is impossible, as the scalar 
perturbations are governed not only by V^(0), but by the flrst derivative of V^(0) as well. 
Knowledge of only the scalar perturbations leaves an undetermined integration constant 
in the non-linear system of reconstruction equations. Another problem is that simple 
potentials usually lead to a nearly exact power-law scalar spectrum with the spectral 
index close to unity. In such a scenario only a very limited amount of information could 
be obtained about high energy physics from astrophysical observations. 

However if a tensor mode can be determined, then one may follow the following 
reconstruction procedure, i?^,, n^., n'^, ... are to be determined from observations. For- 
tunately, parameter estimation from the microwave background has been explored in 
some detail [|T3], n|. We shall use error estimates for Planck assuming polarized detec- 
tors are available, following the analysis of Zaldarriaga et al., [l^. Most analyses have 
assumed that -R* and ri* are the only parameters needed to describe the spectra, 
a recent paper 



m 



In 
Copeland, Grivell, and Liddle have generalized their treatment to 
allow the power spectrum to deviate from scale-invariance. Including extra parameters 
leads to a deterioration in the determination of all the parameters, as it introduces extra 
parameter degeneracies. Fortunately, for most parameters the uncertainty is not much 
increased by including the flrst few derivatives of n |]15[, but the parameter n itself has 
a greatly increased error bar 
An 



0.004 



n, 15 



If a power-law is assumed it can be determined to around 
but including scale dependence increases this error bar by a factor 



of ten or more. Notice that unless one assumes a perfect power-law behavior, this in- 
crease in uncertainty is applicable even if the deviation from power-law behavior cannot 
be detected within the uncertainty. 

From Grivell and Liddle [Hi, an estimate of the relevant uncertainties is 



A(i?,) ~ 0.004 
An ^0.15 



/\{dn/d\nk) ~ 0.04 
^.[d^n/d{\nkf] ~ 0.005 , 



(16) 



Once i?*, y4|, ra*, and n'^ are determined, then it is possible to flnd V{(j)^), V'{ 
and V"{(l)^) from Eq. |15[ Then one can express V(0) as a Taylor series about V^(0*^ 



V^(0) = V^(0,) + V'{ 



,)A0 + ^\/"(0,)A02 



(17) 
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A0 is found from an exact expression connecting changes in with changes in fc 0, 

^^ = ^-^-./^(l + i?.). (18) 

dink An H e-l \ An^ ^ ^ ^ 

Let me illustrate reconstruction by considering two sample potentials. The first po- 
tential, discussed by Lidsey et al, is a power-law potential, V^cj)) = Voexp(— a0/mp;), 
with a ~ 1.6. The potential is shown in the upper left-hand-side of Fig. ^ This potential 
generates a spectral index of n^ = 0.9 and n'^ = 0. It also results in a value of i?* = 0.1. 
So one might guess that precision CBR measurements will determine i?*, n*, and n'^ of 
these values, with uncertainties of Eq. |l^. 

In order to see what sort of reconstructed potential results, one can imagine a model 
universe with R^,, n^,, and tt,'^ generated as Gaussian random variables with mean deter- 
mined by the underlying potential and variance determined by the expected observational 
uncertainties. The result of ten such reconstructions of the potential are shown in the 
upper left-hand figure of Fig. |^. The range of (p is determined by the range of wavenum- 
ber over which one expects to have accurate determinations of the parameters from the 
CBR. Here k^, was chosen to be O.Olh Mpc~^ and the range of k taken to be three decades. 

Also shown in Fig. ^ in the upper right-hand panel is information about the recon- 
struction of the first derivative of the potential. 



The second example of reconstruction was considered in Copeland et al, |16] The 



potential, first considered by Wang et al, ^7j is a potential of the form 



y(0) = AMl--tan-i^) , (19) 

V TT mpij 

which is shown in the lower left-hand side of Fig. ^. 

A fit to a Taylor expansion to the exact spectrum, using k^, as above, yields the 
following results: 

ln/l|(fc,) =3.84 ; i?, = 0.014 ;n, = 0.579 ; 
dn 



dink 



d'^n 
-0.134 



d{lnky 



0.052 . (20) 



From the estimated observational uncertainties of Eq. (0), we see that all these coeffi- 
cients should be successfully determined at high significance and a simple "chi-by-eye" 
demonstrates that the spectrum reconstructed from these data is an adequate fit to the 
observed spectrum. This stresses the point that the more unusual a potential is, the 
more information one is likely to be able to extract about it, though the uncertainties 
on the individual pieces of information may be greater. 

We reconstruct the potential in a region about (p^, = —0.22 (the reconstruction pro- 
gram does not determine 0^,), the width of the region given via Eq. (|THp. The result for 
the reconstructed potential including the effect of observational errors and cosmic vari- 
ance are shown in Fig. ^ where it can immediately be seen that the reconstruction has 
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-'0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.' 
(*-0o)/nip, 




^0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 
(*-0o)/nipi 





Figure 4: Sample reconstructions of two inflationary potentials. The upper two figures 
are the reconstruction of a power-law potential of the form V{(j)) = Voexp(— a0/mpi). 
The light curves are from the true potential, while the heavy curves are ten reconstruc- 
tions. The lower figures are the actual potential and ten reconstructions of a potential 
of the form V{(f)) = A^[l - (2/7r) tan-^(50/mpi)]. 

been very successful in reproducing the main features of the potential while perturbations 
on interesting scales are being developed. 

The uncertainty is dominated by that of y4|,(A;*); although the gravitational waves 
are detectable in this model, it is only about a three-sigma detection and the error bar is 
thus large. Since the overall magnitude of the potential is proportional to A^, the visual 
impression is of a large uncertainty. 

Fortunately, information in combinations of the higher derivatives is more accurately 
determined. Fig. | also shows the reconstruction of V' /V^^'^ with observational errors; 
this combination is chosen as it is independent of the tensors to lowest order. Not only 
is it reconstructed well at the central point, but both the gradient and curvature are 
well fit too, confirming useful information has been obtained about not just V" but V" 
as well, which is only possible because of the extra information contained in the scale- 
dependence of the power spectrum. So rather accurate information is being obtained 
about the potential. 



16 



3 Preheating, Reheating, and Dark Matter 

If the inflaton is completely decoupled, then once inflation ends it will oscillate about the 
minimum of the potential, with the cycle-average of the energy density decreasing as a~^, 
i.e., as a matter-dominated universe. But at the end of inflation the universe is cold and 
frozen in a low-entropy state: the only degree of freedom is the zero-momentum mode 
of the inflaton field. It is necessary to "defrost" the universe and turn it into a "hot" 
high-entropy universe with many degrees of freedom in the radiation. Exactly how this is 
accomplished is still unclear. It probably requires the inflaton field to be coupled to other 
degrees of freedom, and as it oscillates, its energy is converted to radiation either through 
incoherent decay, or through a coherent process involving very complicated dynamics of 
coupled oscillators with time-varying masses. In either case, it is necessary to extract 
the energy from the inflaton and convert it into radiation. 

I will now turn to a discussion of how defrosting might occur. It may be a complicated 
several-step process. I will refer to nonlinear effects in defrosting as "preheating" and 
refer to linear processes as "reheating" |jl8[. 

The possible role of nonlinear dynamics leading to explosive particle production has 
recently received a lot of attention. This process, known as "preheating" 0] may con- 
vert a fair fraction of the inflaton energy density into other degrees of freedom, with 
extremely interesting cosmological effects such as symmetry restoration, baryogenesis, 
or production of dark matter. But the efficiency of preheating is very sensitive to the 
model and the model parameters. 

Perhaps some relic of defrosting, such as symmetry restoration, baryogenesis, or dark 
matter may provide a clue of the exact mechanism, and even shed light on inflation. 

3.1 Defrosting the Universe After Inflation 

3.1.1 Reheating 

In one extreme is the assumption that the vacuum energy of inflation is immediately 
converted to radiation resulting in a reheat temperature Trh- 

A second (and more plausible) scenario is that reheating is not instantaneous, but 
is the result of the slow decay of the inflaton fleld. The simplest way to envision this 
process is if the comoving energy density in the zero mode of the inflaton decays into 
normal particles, which then scatter and thermalize to form a thermal background. It is 
usually assumed that the decay width of this process is the same as the decay width of 
a free inflaton fleld. 

There are two reasons to suspect that the inflaton decay width might be small. The 
requisite flatness of the inflaton potential suggests a weak coupling of the inflaton fleld 
to other flelds since the potential is renormalized by the inflaton coupling to other flelds 
| 19| . However, this restriction may be evaded in supersymmetric theories where the 



nonrenormalization theorem ensures a cancellation between flelds and their superpart- 
ners. A second reason to suspect weak coupling is that in local supersymmetric theories 
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gravitinos are produced during reheating. Unless reheating is delayed, gravitinos will be 
overproduced, leading to a large undesired entropy production when they decay after 



big-bang nucleosynthesis pO |. 



With the above assumptions, the Boltzmann equations describing the redshift and 
interchange in the energy density among the different components is 

/?0 + ?>Hp^ + T^p^ = 

PR + AHpn - T^p^ = , (21) 

where dot denotes time derivative. The dynamics of reheating will be discussed in detail 
in Section B.2.3. 



The reheat temperature is calculated quite easily |T^. After inflation the infiaton field 
executes coherent oscillations about the minimum of the potential. Averaged over several 
oscillations, the coherent oscillation energy density redshifts as matter: p^ oc a~^, where 
a is the Robert son- Walker scale factor. If we denote as pi and aj the total infiaton 
energy density and the scale factor at the initiation of coherent oscillations, then the 
Hubble expansion rate as a function of a is {Mpi is the Planck mass) 



Hia) 



^ft^-l ■ <^^' 



a 



Equating H{a) and T^ leads to an expression for aj/a. Now if we assume that all 
available coherent energy density is instantaneously converted into radiation at this value 
of aj/a, we can define the reheat temperature by setting the coherent energy density. 
Pep = Pi{cii/ciY, equal to the radiation energy density, pr = {ii'^/ 30) g^T^^j, where g^ is 
the effective number of relativistic degrees of freedom at temperature Trh. The result is 

The limit from gravitino overproduction is Trh ^ 10^ to 10^° GeV. 

3.1.2 Preheating 

The main ingredient of the preheating scenario introduced in the early 1990s is the 
nonperturbative resonant transfer of energy to particles induced by the coherently oscil- 
lating infiaton fields. It was realized that this nonperturbative mechanism can be much 
more efficient than the usual perturbative mechanism for certain parameter ranges of 
the theory pl |. 



The basic picture can be seen as follows. Assume there is an infiaton field (p oscillating 
about the minimum of a potential V{(f)) = m'^(f)'^/2. It is convenient to parameterize the 
motion of the infiaton field as 

0(t) = $(t)sinmt. (24) 
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In an expanding universe, even in the absence if interactions the amphtude of the oscil- 
lations of the inflaton field, $, decreases slowly due to the redshift of the momentum. In 
Minkowski space, $ would be a constant in the absence of interactions.. 

Suppose there is a scalar field X with a coupling to the inflaton of g'^cfP'X'^ /2. The 
mode equation for the X field can be written in terms of a redefined variable Xk = Xkd^^"^ 
as 

+ /<l>^(t)sin^(mt) 



a 
Xk + 3-Xfc + 
a 



a2 



Xfe = . (25) 



In Minkowski space a/ a = and $(t) is constant, and the mode equation becomes 
(prime denotes d/dz where z = mt) 

Xl + [Ak - 2qcosi2z)]xkit) = . (26) 

The parameter q depends on the inflaton field oscillation amplitude, and A^ depends on 
the energy of the particle and q: 

^2$2 



Q = 


Am? 




e 


Au -- 


= 2q+ — 



(27) 

When Ak and q are constants, the equation is the Mathieu equation, which exhibits 
resonant mode instabilities for certain values of Ak and q. If Ak and q are constant, then 
there are instability regions where there is explosive particle production. The instability 
regions are shown in Fig. |^ from the paper of Chung [^ . 

Integration of the field equations for the number of particles created in a particular k 
mode is shown in Fig. ^. Explosive growth occurs every time the inflaton passes through 
the origin. 

In an expanding universe, A^ and q will vary in time, but if they vary slowly compared 
to the frequency of oscillations, the effects of resonance will remain. An example of 
particle production where A^ and q vary due to expansion is shown in Fig. |^. Details of 
the parameters and the calculation can be found in Kofman, Linde, and Starobinski p|. 

If the mode occupation number for the X particles is large, the number density per 
mode of the X particles will be proportional to IxfcP- If ^fc and q have the appropriate 
values for resonance, Xk will grow exponentially in time, and hence the number density 
will attain an exponential enhancement above the usual perturbative decay. This period 
of enhanced rate of energy transfer has been called preheating primarily because the 
particles that are produced during this period have yet to achieve thermal equilibrium. 

This resonant amplification leads to an efficient transfer of energy from the inflaton 
to other particles which may have stronger coupling to other particles than the inflaton, 
thereby speeding up the reheating process and leading to a higher reheating temperature 
than in the usual scenario. 

One possible result of preheating is a phase transition caused by the large number of 
soft particles created in preheating |^ . 
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Figure 5: Shaded areas are regions of explosive particle production. The line is A/^ = 2q 
(recall A^ = 2q + k"^ /m?). This figure is from Chung [^ . 




Figure 6: Explosive particle production if Ak and q are constant. This figure is from 
Kofman, Linde, and Starobinski HI. 
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Inn^, 




Figure 7: Explosive particle production if Ak and q are not constant. The lower 
figure shows integration for a longer time. These figures are from Kofman, Linde, and 
Starobinski M. 
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Another interesting feature is that particles of mass larger than the infiaton mass can 
be produced through this coherent resonant effect. This has been exploited to construct 



a baryogenesis scenario |24] in which the baryon number violating bosons with masses 



larger than the infiaton mass are created through the resonance mechanism. 

3.2 Dark Matter 

There are many reasons to believe the present mass density of the universe is dominated 
by a weakly interacting massive particle (wiMP), a fossil relic of the early universe. The- 
oretical ideas and experimental efforts have focused mostly on production and detection 
of thermal relics, with mass typically in the range a few GeV to a hundred GeV. Here, 
I will review scenarios for production of nonthermal dark matter. Since the masses of 
the nonthermal WIMPS are in the range 10^^ to 10^^ GeV, much larger than the mass 
of thermal wimpy WIMPS, they may be referred to as wimpzillas. In searches for dark 
matter it may be well to remember that "size does matter." 

3.2.1 Thermal Relics — Wimpy WIMPS 

It is usually assumed that the dark matter consists of a species of a new, yet undiscovered, 
massive particle, traditionally denoted by X. It is also often assumed that the dark 
matter is a thermal relic, i.e., it was in chemical equilibrium in the early universe. 

A thermal relic is assumed to be in local thermodynamic equilibrium (lte) at early 
times. The equilibrium abundance of a particle, say relative to the entropy density, 
depends upon the ratio of the mass of the particle to the temperature. Define the 
variable Y = nx/s, where nx is the number density of WIMP X with mass Mx, and 
s ~ T^ is the entropy density. The equilibrium value of Y, Yeq, is proportional to 
exp(— x) for X ^ 1, while Yeq ~ constant for x <^ 1, where x = Mx/T. 

A particle will track its equilibrium abundance as long as reactions which keep the 
particle in chemical equilibrium can proceed rapidly enough. Here, rapidly enough means 
on a timescale more rapid than the expansion rate of the universe, H. When the reaction 
rate becomes smaller than the expansion rate, then the particle can no longer track its 
equilibrium value, and thereafter Y is constant. When this occurs the particle is said to 
be "frozen out." A schematic illustration of this is given in Fig. p. 

The more strongly interacting the particle, the longer it stays in lte, and the smaller 
its eventual freeze-out abundance. Conversely, the more weakly interacting the particle, 
the larger its present abundance. The freeze-out value of Y is related to the mass of the 
particle and its annihilation cross section (here characterized by ao) by \TE 



Y oc — . (28) 

Mxmpiao 

Since the contribution to Q is proportional to MxUx, which in turn is proportional to 
MxY, the present contribution to Q from a thermal relic roughly is independent of its 
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10' 108 

X = Mx/T 



Figure 8: A thermal relic starts in LTE at T ^ Mx- When the rates keeping the relic 
in chemical equilibrium become smaller than the expansion rate, the density of the relic 
relative to the entropy density freezes out. 



mass,[] and depends only upon the annihilation cross section. The cross section that 
results in Vtxh'^ ~ 1 is of order 10~^^cm^, i.e., of the order of the weak scale. This is one 
of the attractions of thermal relics. The scale of the annihilation cross section is related 
to a known mass scale. 

The simple assumption that dark matter is a thermal relic is surprisingly restric- 
tive. The largest possible annihilation cross section is roughly M^^. This implies that 
large-mass wimps would have such a small annihilation cross section that their present 
abundance would be too large. Thus one expects a maximum mass for a thermal WIMP, 
which turns out to be a few hundred TeV [^ 



The standard lore is that the hunt for dark matter should concentrate on particles 
with mass of the order of the weak scale and with interaction with ordinary matter on 
the scale of the weak force. This has been the driving force behind the vast effort in dark 
matter direct detection. 

In view of the unitarity argument, in order to consider thermal wimpzillas, one 
must invoke, for example, late-time entropy production to dilute the abundance of these 
supermassive particles |^|, rendering the scenario unattractive. 



3.2.2 Nonthermal Relics— WIMPZILLAS 

There are two necessary conditions for the wimpzilla scenario. First, the wimpzilla 
must be stable, or at least have a lifetime much greater than the age of the universe. 



''To first approximation the relic dependence depends upon the mass only indirectly through the 
dependence of the annihilation cross section on the mass. 
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This may result from, for instance, supersymmetric theories where the breaking of su- 
persymmetry is communicated to ordinary sparticles via the usual gauge forces [^. In 



particular, the secluded and the messenger sectors often have accidental symmetries anal- 
ogous to baryon number. This means that the lightest particle in those sectors might be 



stable and very massive if supersymmetry is broken at a large scale [53]. Other natural 



candidates arise in theories with discrete gauge symmetries ||2^ and in string theory and 
M theory m,m. 



It is useful here to note that wimpzilla decay might be able to account for ultra-high 
energy cosmic rays above the Greisen-Zatzepin-Kuzmin cutoff [RSI, RBI. A wimpy little 



thermal relic would be too light to do the job, a wimpzilla is needed. 

The second condition for a wimpzilla is that it must not have been in equilibrium 
when it froze out {i.e., it is not a thermal relic), otherwise Vtxh? would be much larger 
than one. A sufficient condition for nonequilibrium is that the annihilation rate (per 
particle) must be smaller than the expansion rate: nxcif | < H, where a\v\ is the 
annihilation rate times the M0ller flux factor, and H is the expansion rate. Conversely, 
if the dark matter was created at some temperature T^, and Vtxh? < 1, then it is easy to 
show that it could not have attained equilibrium. To see this, assume X's were created 
in a radiation-dominated universe at temperature T^. Then Vtxh? is given by 

^xh^ = n^h\TjTo)mxnxiT,)/p^iT,) , (29) 

where Tq is the present temperature. Using the fact that p-yiT^) = H{T^)MpiT^, 
nxiT^)/H{T^) = (nx/^'y)ToMpiT^/Mx- One may safely take the limit cr|i;| < Mx'^, 
so nxiTt,)a\v\/H{T^:) must be less than {fix/^'y)ToMpiT^,/Mx- Thus, the requirement 
for nonequilibrium is 

This implies that if a nonrelativistic particle with Mx ~ 200 TeV was created at T* < Mx 
with a density low enough to result in Qx ~ 1, then its abundance must have been so small 
that it never attained equilibrium. Therefore, if there is some way to create wimpzillas 
in the correct abundance to give fix ~ 1, nonequilibrium is automatic. Examples of 
WIMPZILLA evolution and freezeout are shown in Fig. |^. 

Any WIMPZILLA production scenario must meet these two criteria. Before turning to 
several wimpzilla production scenarios, it is useful to estimate the fraction of the total 
energy density of the universe in wimpzillas at the time of their production that will 
eventually result in fi ~ 1 today. 

The most likely time for wimpzilla production is just after inflation. The first step 
in estimating the fraction of the energy density in wimpzillas is to estimate the total 
energy density when the universe is "reheated" after infiation. 

Consider the calculation of the reheat temperature, denoted as TpH. The reheat 
temperature is calculated by assuming an instantaneous conversion of the energy density 
in the infiaton field into radiation when the decay width of the infiaton energy, F^, is 
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Mx/T 



Mx/T 



Figure 9: A nonthermal relic typically has abundance much less than lte at T ^ Mx- 
Here are two examples of the evolution of a nonthermal relic. In the left panel Y -C Yeq 
at freezeout, while in the right panel Y ^ Yeq at freezeout. Again, the solid curve is 
the LTE abundance while the dashed curve is the actual abundance. 

equal to if, the expansion rate of the universe. The derivation of the reheat temperature 
was given in Section p.l.l| . 

Now consider the WIMPZILLA density at reheating. Suppose the wimpzilla never 
attained lte and was nonrelativistic at the time of production. The usual quantity Vtxh^ 
associated with the dark matter density today can be related to the dark matter density 
when it was produced. First write 



Px{to) Px{tRH) (Trh 



pRito) pRitRu) V Tc 



(31) 



where pr denotes the energy density in radiation, px denotes the energy density in 
the dark matter, Tj^j is the reheat temperature, Tq is the temperature today, to de- 
notes the time today, and t^H denotes the approximate time of reheating.^ To obtain 
Px(tRH)/PR(tRH), one must determine when X particles are produced with respect to 
the completion of reheating and the effective equation of state between X production 
and the completion of reheating. 

At the end of inflation the universe may have a brief period of matter domination 
resulting either from the coherent oscillations phase of the infiaton condensate or from 
the preheating phase PT|. If the X particles are produced at time t = t^ when the 
de Sitter phase ends and the coherent oscillation period just begins, then both the X 
particle energy density and the infiaton energy density will redshift at approximately 
the same rate until reheating is completed and radiation domination begins. Hence, the 



^More specifically, this is approximately the time at which the universe becomes radiation dominated 
after inflation. 
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ratio of energy densities preserved in this way until the time of radiation domination is 

Px{tRH) Svr px{te) 



PR{tRH) 3 MliH^{t 



e/ 



(32) 



where Mpi ^ 10^^ GeV is the Planck mass and most of the energy density in the universe 
just before time tun is presumed to turn into radiation. Thus, using Eq. ^, one may 
obtain an expression for the quantity fix = Px(to)/pc(to), where pc(^o) = SHqMpi/Sti 
and Hq = 100 h km sec~^ Mpc~^: 

"" \To J 3 \MpJ MpiH\Q ^ ' 

Here flph'^ ~ 4.31 x 10~^ is the fraction of critical energy density in radiation today and 
nx is the density of X particles at the time when they were produced. 

Note that because the reheating temperature must be much greater than the temper- 
ature today {Trh/Tq ^ 4.2 x 10^"^), in order to satisfy the cosmological bound Qxh"^ ~ 1, 
the fraction of total wimpzilla energy density at the time when they were produced must 
be extremely small. One sees from Eq. |3|that ilxh'^ ~ IO^^Trh /lO^GeV){pxite) / pite))- 
It is indeed a very small fraction of the total energy density extracted in wimpzillas. 

This means that if the wimpzilla is extremely massive, the challenge lies in creating 



very few of them. Gravitational production discussed in Section |3.2.3| naturally gives 
the needed suppression. Note that if reheating occurs abruptly at the end of inflation, 
then the matter domination phase may be negligibly short and the radiation domination 
phase may follow immediately after the end of inflation. However, this does not change 
Eq. ||. 

3.2.3 Wimpzilla Production 

3. 2. 3. a Gravitational Production 
First consider the possibility that wimpzillas are produced in the transition between 
an inflationary and a matter-dominated (or radiation-dominated) universe due to the 
"nonadiabatic" expansion of the background spacetime acting on the vacuum quantum 
fluctuations 133 . 



The distinguishing feature of this mechanism is the capability of generating particles 
with mass of the order of the inflaton mass (usually much larger than the reheating 
temperature) even when the particles only interact extremely weakly (or not at all) with 
other particles and do not couple to the inflaton. They may still be produced in sufficient 
abundance to achieve critical density today due to the classical gravitational effect on 
the vacuum state at the end of inflation. More specifically, if 0.04 < Mx/Hj < 2, where 
Hj ~ m^ ~ 10^'^GeV is the Hubble constant at the end of inflation [m^ is the mass 
of the inflaton), wimpzillas produced gravitationally can have a density today of the 
order of the critical density. This result is quite robust with respect to the "flne" details 
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Figure 10: The contribution of gravitationally produced wimpzillas to Qxh"^ as a 
function of Mx/Hj. The shaded area is where thermahzation may occur if the annihila- 
tion cross section is its maximum value. Also shown is the contribution assuming that 
the WIMPZILLA is present at the end of inflation with a temperature T = ilf//27r. 



of the transition between the inflationary phase and the matter-dominated phase, and 
independent of the coupling of the wimpzilla to any other particle. 

Conceptually, gravitational wimpzilla production is similar to the inflationary gen- 
eration of gravitational perturbations that seed the formation of large scale structures. 
In the usual scenarios, however, the quantum generation of energy density fluctuations 
from inflation is associated with the inflaton field that dominated the mass density of 
the universe, and not a generic, sub-dominant scalar field. Another difference is that 
the usual density fluctuations become larger than the Hubble radius, while most of the 
WIMPZILLA perturbations remain smaller than the Hubble radius. 

There are various inequivalent ways of calculating the particle production due to in- 
teraction of a classical gravitational field with the vacuum (see for example ||3^, ^ W^)- 
Here, I use the method of finding the Bogoliubov coefficient for the transformation be- 
tween positive frequency modes defined at two different times. For Mx/Hj < 1 the 
results are quite insensitive to the differentiability or the fine details of the time depen- 
dence of the scale factor. For 0.03 < Mx/Hj < 10, all the dark matter needed for closure 
of the universe can be made gravitationally, quite independently of the details of the 
transition between the infiationary phase and the matter dominated phase. 

Start with the canonical quantization of the X field in an action of the form (with 
metric ds'^ = dt^ — a'^{t)dy? = a?{ri) [drf — dy?] where 7] is conformal time) 



S 



dt d^x— X^ 



{vxy 



MJ^X^ 



iRX' 



(34) 
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where R is the Ricci scalar. After transforming to confornial time coordinate, use the 
mode expansion 



XM 



d\ 



(27r)3/2a(r7) 



akh{r])e'^-^ + alhl{7])e-'^-^ , (35) 



where because the creation and annihilation operators obey the commutator [a^^, a^ = 
5''^'*(ki — k2), the hkS obey a normalization condition hkhj^ — h'f^h\ = i to satisfy the 
canonical field commutators (henceforth, all primes on functions of rj refer to derivatives 
with respect to rj). The resulting mode equation is 

Kiv) + wlir^)h,{r^) = 0, (36) 

where 

wl = k'^ + Mla^ + (6^ - l)a"/a . (37) 

The parameter ^ is 1/6 for conformal coupling and for minimal coupling. From now 
on, ^ = 1/6 for simplicity but without much loss of generality. By a change in variable 
7] -^ k/a, one can rewrite the differential equation such that it depends only on H{ri), 
H'{ri)/k, k/a{ri), and Mx- Hence, the parameters Hj and aj correspond to the Hubble 
parameter and the scale factor evaluated at an arbitrary conformal time rjj, which can 
be taken to be the approximate time at which Xs are produced {i.e., rji is the conformal 
time at the end of inflation). 



One may then rewrite Eq. |36| as 



Ml 



h'm+[k' + ^~a']h~M = 0, (38) 



where fj = rjajHi, a = a/ai, and k = k/{aiHj). For simplicity of notation, drop all 
the tildes. This differential equation can be solved once the boundary conditions are 
supplied. 

The number density of the wimpzillas is found by a Bogoliubov transformation 
from the vacuum mode solution with the boundary condition at i] = tjq (the initial time 
at which the vacuum of the universe is determined) into the one with the boundary 
condition at r^ = r^i (any later time at which the particles are no longer being created). 
Tjo will be taken to be — cxd while rji will be taken to be at +cx3. Defining the Bogoliubov 
transformation as hP^ij]) = akhJj^ijj) + Pkhl^^iv) (the superscripts denote where the 
boundary condition is set), the energy density of produced particles is 

/ 1 \^ rc<y Hk ~ 
pxivi) = Mxnxim) = MxHf i^-^^j J^ j^.k^W^ (39) 

where one should note that the number operator is defined at rji while the quantum 
state (approximated to be the vacuum state) defined at tjq does not change in time in 
the Heisenberg representation. 
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Figure 11: The evolution of the Bogohubov coefficient with conformal time for several 
wavenumbers. rj = rjj corresponds to the end of the inflationary era. 



As one can see from Eq. |3^, the input parameter is Mx/Hj. One must also spec- 
ify the behavior of a{ri) near the end of inflation. In Fig. |10| (from |3|]), I show the 
resulting values of Vtxh? as a function of Mx/Hj assuming the evolution of the scale 
factor smoothly interpolates between exponential expansion during inflation and either 
a matter-dominated universe or radiation-dominated universe. The peak at Mx/Hj ~ 1 
is similar to the case presented in Ref. ||38|. As expected, for large Mx/Hj, the number 
density falls off faster than any inverse power of Mx/Hj. 

Now most of the action occurs around the transition from inflation to the matter- 



dominated or radiation-dominated universe. This is shown in Fig. [T^. Also from Fig. Ill 
one can see that most of the particles are created with wavenumber of order Hi. 

To conclude, there is a significant mass range (0.03iJ/ to lOHj, where Hj ~ lO^^GeV) 
for which WIMPZILLAS will have critical density today regardless of the fine details of 
the transition out of inflation. Because this production mechanism is inherent in the 
dynamics between the classical gravitational field and a quantum field, it needs no fine 
tuning of field couplings or any coupling to the inflaton field. However, only if the 
particles are stable (or sufficiently long lived) will these particles give contribution of the 
order of critical density. 

3.2.3.b Production During Reheating 
Another attractive origin for WIMPZILLAS is during the defrosting phase after infla- 
tion. It is important to recall that it is not necessary to convert a significant fraction of 
the available energy into massive particles; in fact, it must be an infinitesimal amount. 
I now will discuss how particles of mass much greater than T^^h may be created in the 
correct amount after inflation in reheating 
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In one extreme is the assumption that the vacuum energy of inflation is immediately 
converted to radiation resulting in a reheat temperature Trh. In this case Qx can 
be calculated by integrating the Boltzmann equation with initial condition Nx = at 
T = Trh- One expects the X density to be suppressed by exp(— 2Mx/rR//); indeed, one 
finds VLx ~ 1 for Mx/Trh ~ 25 + 0.5 ln(m^((T|t'|)), in agreement with previous estimates 
3^ that for Trh ~ lO^GeV, the wimpzilla mass would be about 2.5 x 10^°GeV. 

It is simple to calculate the wimpzilla abundance in the slow reheating scenario. It 
will be important to keep in mind that what is commonly called the reheat temperature, 
Trh, is not the maximum temperature obtained after inflation. The maximum temper- 
ature is, in fact, much larger than Trh- The reheat temperature is best regarded as 
the temperature below which the universe expands as a radiation-dominated universe, 
with the scale factor decreasing as g* T~^. In this regard it has a limited meaning 
40]. One implication of this is that it is incorrect to assume that the maximum 



abundance of a massive particle species produced after inflation is suppressed by a factor 

ofexp(-M/THH). 

To estimate WIMPZILLA production in reheating, consider a model universe with 
three components: inflaton field energy, p^, radiation energy density, Pr, and WIMPZILLA 
energy density, px- Assume that the decay rate of the inflaton field energy density is F,^. 
Also assume the wimpzilla lifetime is longer than any timescale in the problem (in fact 
it must be longer than the present age of the universe). Finally, assume that the light 
degrees of freedom are in local thermodynamic equilibrium. 

With the above assumptions, the Boltzmann equations describing the redshift and 
interchange in the energy density among the different components is [cf. Eq. (PT|)] 



P0 + 3if p0 + F 
PR + AH PR - F 





Mx 



px + "iHpx 
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2 ( EQ\ 
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[pI^) 







(40) 



where dot denotes time derivative. As already mentioned, (o"|t'|) is the thermal average 
of the X annihilation cross section times the M0ller flux factor. The equilibrium energy 
density for the X particles, Px ■, is determined by the radiation temperature, T = 
(30pK/vrV)i/l 

It is useful to introduce two dimensionless constants, a^ and ax-, defined in terms of 
F<^ and (cr|f |) as 

F,^ = a^M^ (,a\v\) = axM^^ . (41) 

For a reheat temperature much smaller than M^, F^ must be small. From Eq. (^), the 
reheat temperature in terms of ax and Mx is Trh ~ a^ J M^Mpi. For M^ = lO^'^'^GeV, 
a^ must be smaller than of order 10~^^. On the other hand, ax may be as large as of 
order unity, or it may be small also. 
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It is also convenient to work with dimensionless quantities that can absorb the effect 
of the expansion of the universe. This may be accomphshed with the definitions 

$ = P4,M^'a^ ; R = PRa^ ; X = pxM^'a^ . (42) 

It is also convenient to use the scale factor, rather than time, for the independent variable, 
so one may define a variable x = aM^. With this choice the system of equations can be 
written as (prime denotes d/dx) 

$' = -ci , "^ $ 

yJ^X + R 



1 —1 

R' = ci ^^^ $ + C2 ^^^ f X^ - X^ 



^J^x + R ^J^x + R 



EQ 



X '^ 



X' = -C3 / (X^ - Xlo) . (43) 

^^x + R \ ^QJ ^ ' 

The constants Ci, C2, and C3 are given by 

^^ = V^^"^ '' = ''J^^ '' = ''Wx- ^''^ 

Xeq is the equilibrium value of X, given in terms of the temperature T as (assuming a 
single degree of freedom for the X species) 

The temperature depends upon R and g^,, the effective number of degrees of freedom in 
the radiation: 

X.fiLV"^^ (46) 

It is straightforward to solve the system of equations in Eq. ( ^31) with initial conditions 
a.t X = xj of R{xj) = X{xj) = and ^{xj) = $7. It is convenient to express p^{x = xj) 
in terms of the expansion rate at x/, which leads to 

MMlMl ' ^ ' 

The numerical value of x/ is irrelevant. 

Before solving numerically the system of equations, it is useful to consider the early- 
time solution for R. Here, early times means H ^ T^, i.e., before a significant fraction of 
the comoving coherent energy density is converted to radiation. At early times $ ~ $/, 
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Figure 12: The evolution of energy densities and T jMx as a function of the scale factor. 
Also shown is X/ Xeq- 

and i? ~ X ~ 0, so the equation for R' becomes R' = ciX^/^$/ . Thus, the early time 
solution for R is simple to obtain: 



2 / 
R ~ -ci (x 



5/2 



xf) $f 



fi^>r, 



(48) 



Now express T in terms of R to yield the early-time solution for T: 



T 

Ml 



12 



1/4 



c^ 



i,' 



X J 



X \ ^/^ 

Xj , 



X 
Xi 



1/4 



[i^>r„ 



Thus, T has a maximum value of 



(49) 
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(50) 



which is obtained at x/xi = (8/3)^/^ = 1.48. It is also possible to express a^ in terms of 
TjiH and obtain 



Tmax _ r.yy I 9 



-RH 



hn^g* 



HiMpi 

^ RH 



1/4 



(51) 



For an illustration, in the simplest model of chaotic inflation Hj ~ M^ with M^ 
lO^^GeV, which leads to Tmax/Trh ~ 10=^(200/^,)^/^ for Trh = lO^GeV. 
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We can see from Eq. p8|) that for x/xi > 1, in the early-time regime T scales as 
a~^/®, which implies that entropy is created in the early-time regime [0. So if one is 
producing a massive particle during reheating it is necessary to take into account the 
fact that the maximum temperature is greater than Trh-, and that during the early-time 
evolution, T oc a~^/^. 

An example of a numerical evaluation of the complete system in Eq. (|43|) is shown in 
Fig. |I2] (from [§^). The model parameters chosen were M^ = lO^^GeV, a^ = 2 x 10~^^, 
Mx = 1.15 X lO^^GeV, ax = 10^^, and (?* = 200. The expansion rate at the beginning 
of the coherent oscillation period was chosen to be Hj = M^. These parameters result 
in Trh = lO^GeV and ^xh"^ = 0.3. 



Figure |I2| serves to illustrate several aspects of the problem. Just as expected, the 
comoving energy density of {i.e., a?p^) remains roughly constant until F,^ ~ if, which 
for the chosen model parameters occurs around a/ai ^ 5 x 10^. But of course, that does 
not mean that the temperature is zero. Notice that the temperature peaks well before 
"reheating." The maximum temperature, Tmax = lO^^GeV, is reached at a/aj slightly 
larger than unity (in fact at a/ai = 1.48 as expected), while the reheat temperature, 
Trh = lO^GeV, occurs much later, around a/aj ~ 10*. Note that Tmax — IO'^Trj^ in 
agreement with Eq. (0). 

From the figure it is clear that X <ti Xeq at the epoch of freeze out of the comoving 
X number density, which occurs around a/aj ~ 10^. The rapid rise of the ratio after 
freeze out is simply a reflection of the fact that X is constant while Xeq decreases 
exponentially. 

A close examination of the behavior of T shows that after the sharp initial rise of the 
temperature, the temperature decreases as a~'^/^ [as follows from Eq. (^91)] until H ~ F,^, 
and thereafter T oc a~^ as expected for the radiation-dominated era. 

For the choices of M^, a^, g^., and ax used for the model illustrated in Fig. |1^, 
Qxh"^ = 0.3 for Mx = 1-15 x lO^^GeV, in excellent agreement with the mass predicted 



by using an analytic estimate for the result |39 



...-MiH.|)(|-)-^'^(^)'. 

Here again, the results have also important implications for the conjecture that ultra- 
high cosmic rays, above the Greisen-Zatsepin-Kuzmin cut-off of the cosmic ray spectrum, 
may be produced in decays of superheavy long-living particles ||3l|, |3^, |3^, ^. In order to 



produce cosmic rays of energies larger than about 10^'^ GeV, the mass of the X-particles 
must be very large, Mx ~ 10^^ GeV and their lifetime tx cannot be smaller than the 
age of the Universe, tx ~ 10^° yr. With the smallest value of the lifetime, the observed 
flux of ultra-high energy cosmic rays will be reproduced with a rather low density of 
X-particles, fix ~ 10~^^. It has been suggested that X-particles can be produced in the 
right amount by usual collisions and decay processes taking place during the reheating 
stage after inflation if the reheat temperature never exceeded Mx [0- Again, assuming 
naively that that the maximum number density of a massive particle species X produced 

33 



after inflation is suppressed by a factor of {Mx /TrhY^"^ Q'^v{~Mx /Trh) with respect to 
the photon number density, one concludes that the reheat temperature T^u should be 
in the range 10^^ to lO^^GeV [32]. This is a rather high value and leads to the gravitino 
problem in generic supersymmetric models. This is one reason alternative production 
mechanisms of these superheavy X-particles have been proposed [^ ^ ^ . However, 



our analysis show that the situation is much more promising. Making use of Eq. (p2D, 
the right amount of X-particles to explain the observed ultra-high energy cosmic rays is 
produced for 

IeL\ ^ (±X''' ( ^^ ^ f53) 

IQiOGeVy V200y UQiSGeVy' ^ ' 

where it has been assumed that (o"|f |) ~ Mx'^. Therefore, particles as massive as 10^^ 
GeV may be generated during the reheating stage in abundances large enough to explain 
the ultra-high energy cosmic rays even if the reheat temperature satisfies the gravitino 
bound. 

3.2.3.C Production During Preheating 

Another way to produce WIMPZILLAS after inflation is in a preliminary stage of 
reheating called "preheating" ||2T|, where nonlinear quantum effects may lead to an 
extremely effective dissipational dynamics and explosive particle production. 

As discussed in Section p.l.2| , particles can be created in a broad parametric resonance 
with a fraction of the energy stored in the form of coherent inflaton oscillations at the 
end of inflation released after only a dozen oscillation periods. A crucial observation 
for our discussion is that particles with mass up to 10^^ GeV may be created during 
preheating ||4^, ^, ^, and that their distribution is nonthermal. If these particles are 



stable, they may be good candidates for wimpzillas [|22 |. 

Interestingly enough, what was found [^ is that in the context of a slow-roll inflation 
with the potential V{(j)) = m^0^/2 with the inflaton coupling oi g'^(f)'^X'^/2, the resonance 
phenomenon is mostly irrelevant to WIMPZILLA production because too many particles 
would be produced if the resonance is effective. For the tiny amount of energy conversion 
needed for wimpzilla production, the coupling g^ must be small enough (for a fixed Mx) 
such that the motion of the inflaton field at the transition out of the inflationary phase 
generates just enough nonadiabaticity in the mode frequency to produce wimpzillas. 
The rest of the oscillations, damped by the expansion of the universe, will not contribute 
significantly to wimpzilla production as in the resonant case. In other words, the quasi- 
periodicity necessary for a true resonance phenomenon is not present in the case when 
only an extremely tiny fraction of the energy density is converted into wimpzillas. Of 
course, if the energy scales are lowered such that a fair fraction of the energy density can 
be converted to wimpzillas without overdosing the universe, this argument may not 
apply. 



The main finding of a detailed treatment |22| is that wimpzillas with a mass as large 



as 10 if/, where Hj is the value of the Hubble expansion rate at the end of inflation, can 
be produced in sufficient abundance to be cosmologically significant today. 
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Figure 13: A graph of Vlxh"^ versus Mx/Hj for gMpi/Hj = 10^. The sohd curve is a 



numerical result, while the dashed and dotted curves are analytic approximations [22 



If the WIMPZILLA is coupled to the inflaton 
equation in Eq. BTI is now changed to 



by a term g (f) X /2, then the mode 



ujl + e+{M], + g'<f)' 



(54) 



again taking C, = 1/6. 

The procedure to calculate the wimpzilla density is the same as in Section |3.2.3 . 
Now, in addition to the parameter Mx/Hj, there is another parameter gMpi/Hj. Now 
in large-field models Hj ~ lO^^GeV, so Mpi/Hj might be as large as 10^. The choice of 
g = 10-3 would yield gMpi/Hi = 10^. 

Fig. |l^ (from p2[) shows the dependence of the WIMPZILLA density upon Mx/Hi 
for the particular choice gMpi/Hj = 10^. This would correspond to ^f ~ 1 in large- 
field inflation models where Mpi/Hj = 10^, about the largest possible value. Note that 
Qx ~ 1 obtains for Mx/Hj ^ 10^. The dashed and dotted curves are two analytic 
approximations discussed in E^, while the solid curve is the numerical result. The 
approximations are in very good agreement with the numerical results. 

Fig. |T^ (from |2^) shows the dependence of the wimpzilla density upon gMpi/Hj. 
For this graph Mx/Hj was chosen to be unity. This figure illustrates the fact that the 
dependence of Vtxh? on gMpi/Hj is not monotonic. For a detailed explanation of this 



curious effect, see the paper of Chung |22 



3.2.3.d Production in Bubble Collisions 
WIMPZILLAS may also be produced |^| if inflation is completed by a first-order phase 



transition [^ , in which the universe exits from a false- vacuum state by bubble nucleation 
W^ . When bubbles of true vacuum form, the energy of the false vacuum is entirely 
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Figure 14: An illustration of the nonmonotonic behavior of the particle density produced 
with the variation of the coupling constant. The value of Mx/Hj is set to unity p^ . 



transformed into potential energy in the bubble walls. As the bubbles expand, more and 
more of their energy becomes kinetic as the walls become highly relativistic. 

In bubble collisions the walls oscillate through each other [^ and their kinetic energy 



is dispersed into low-energy scalar waves |^, |4^. We are interested in the potential 
energy of the walls, Mp = AirrjR'^, where r] is the energy per unit area of a bubble wall 
of radius R. The bubble walls can be visualized as a coherent state of inflaton particles, 
so the typical energy E of the products of their decays is simply the inverse thickness of 
the wall, E ~ A~^. If the bubble walls are highly relativistic when they collide, there 
is the possibility of quantum production of nonthermal particles with mass well above 
the mass of the inflaton field, up to energy A^^ = 7M^, with 7 the relativistic Lorentz 
factor. 

Suppose for illustration that the wimpzilla is a fermion coupled to the inflaton field 
by a Yukawa coupling g(j)XX. One can treat (the bubbles or walls) as a classical, 
external field and the wimpzilla as a quantum field in the presence of this source. 
The number of WIMPZILLAS created in the collisions from the wall potential energy is 
Nx ~ fxMp/Mx, where fx parametrizes the fraction of the primary decay products in 
WIMPZILLAS. The fraction fx will depend in general on the masses and the couplings of 
a particular theory in question. For the Yukawa coupling g, it is fx — 5'^ln {'^M^/2Mx) 
p8| , ^ . WIMPZILLAS may be produced in bubble collisions out of equilibrium and never 
attain chemical equilibrium. Even with Trh as low as 100 GeV, the present wimpzilla 
abundance would be Vtx ~ 1 if (7 ~ 10~^a^/^. Here a~^ <^ 1 is the fraction of the bubble 
energy at nucleation in the form of potential energy at the time of collision. This simple 
analysis indicates that the correct magnitude for the abundance of WIMPZILLAS may be 
naturally obtained in the process of reheating in theories where inflation is terminated 
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by bubble nucleation. 

3.2.4 Wimpzilla Conclusions 

In this talk I have pointed out several ways to generate nonthermal dark matter. All of 
the methods can result in dark matter much more massive than the feeble little weak- 
scale mass thermal relics. The nonthermal dark matter may be as massive as the GUT 
scale, truly in the wimpzilla range. 

The mass scale of the wimpzillas is determined by the mass scale of inflation, more 
exactly, the expansion rate of the universe at the end of inflation. For large-field inflation 
models, that mass scale is of order lO^^GeV. For small-field inflation models, it may be 
less, perhaps much less. 

The mass scale of inflation may one day be measured! In addition to scalar density 
perturbations, tensor perturbations are produced in inflation. The tensor perturbations 
are directly proportional to the expansion rate during inflation, so determination of a 
tensor contribution to cosmic background radiation temperature fluctuations would give 
the value of the expansion rate of the universe during inflation and set the scale for the 
mass of the wimpzilla. 

Undoubtedly, other methods for wimpzilla production will be developed. But per- 
haps even with the present scenarios one should start to investigate methods for wim- 
pzilla detection. While wimpy wimps must be color singlets and electrically neutral, 
wimpzillas may be endowed with color and electric charge. This should open new 
avenues for detection and exclusion of wimpzillas. The lesson here is that wimpzillas 
may surprise and be the dark matter, and we may learn that size does matter! 
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